Quantitative test of mean-field description of a trapped two-dimensional Bose gas 
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We investigate the accuracy of two mean-field theories of the trapped two-dimensional Bose gas at predicting 
transition region properties by comparison to non-perturbative classical field calculations. To make these com- 
parisons we examine the density profiles and the predictions for the Berezinskii-Kosterlitz-Thouless superfluid 
transition temperature over a parameter range in which the degree of thermal activation in the tightly trapped 
direction varies considerably. These results present an important test of these mean-field theories, and provide a 
characterization of their typical accuracy. 
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I. INTRODUCTION 

Evidence for the Berezinskii-Kosterlitz-Thouless (BKT) 
superfluid transition in a trapped two-dimensional (2D) Bose 
gas has been reported in a series of experiments by the ENS, 
NIST and Harvard groups |H0,[l,0,|5t]. Comprehensive theo- 
retical descriptions of the trapped system have been provided 
by classical field (c-field) and quantum Monte Carlo (QMC) 
methods Ja 0, H, Hfl, building on the earlier theoretical work 
for the homogeneous Bose gas 

It is desirable to have a simple theoretical description for 
the 2D trapped Bose gas that can be more easily employed 
than c-field or QMC methods. In 2008, Holzmann, Chevallier 
and Krauth (HCK) 111 311 presented an important step towards 
this goal: a semiclassical mean-field theory they used to pre- 
dict the BKT superfluid transition temperature, Tb kt- This 
work differed from earlier mean-field work (e.g. see 
which aimed to understand the low temperature properties of 
the gas, and stability of the Bose condensed phase against 
thermal fluctuations. The HCK theory instead was applied to 
the high temperature regime (without condensate) to calculate 
the system density profile. From this the BKT superfluid tran- 
sition could be identified by locating the temperature at which 
the peak phase space density of the gas obtained the critical 
value determined from studies of the uniform system ifloll . In 
lfl3ll several comparisons were made between the HCK the- 
ory and full QMC calculations showing good agreement for 
the density profiles and a prediction of Tbkt- 

In previous work lfl6ll we developed a more complete semi- 
classical mean-field theory and a consistent procedure identi- 
fying Tbkt- For clarity we refer to this as the full semiclassi- 
cal (FSC) mean-field theory. While the FSC theory reduces to 
the HCK with some ad hoc simplifications, in the degenerate 
regime the predictions of the two theories differ appreciably. 
Indeed, the direct comparisons between the two mean-field 
theories made in [ 16] demonstrated that FSC always predicts 
lower values for both the system density profile and Tbkt- 
Notably, for large systems the differences between the predic- 
tions was significant. It was argued that the good agreement 
between the HCK mean-field theory and the QMC calcula- 
tions was likely fortuitous. 

Due to the increase in interest in the 2D trapped Bose gas 
it is necessary to better understand the reliability of these 



mean-field theories for making accurate predictions. In this 
work we provide quantitative tests of both mean-field theo- 
ries against c-field calculations. We present results spanning 
a broad regime in which the system crosses from being pure- 
2D (negligible thermal activation in the tightly trapped direc- 
tion) to quasi-2D (with appreciable thermal activation). We 
compare the theories through their predictions for the density 
profiles and the superfluid transition temperature. Our results 
provide a benchmark test for the use of mean-field theory, and 
caution against its use as a quantitative tool. Indeed, our re- 
sults show that the good agreement observed between HCK 
mean-field theory and QMC in Ref. 111311 was due to the quasi- 
2D regime considered there, and such good agreement is not 
found in the pure-2D regime. 



II. QUASI-2D MEAN-FIELD THEORY 

In this section we review the mean-field theory devel- 
oped in Refs. (H, [T^]. We consider a system of ultra-cold 
bosons confined in a harmonic potential with trap frequencies 
{u> x ,u>y,uj z }, interacting via a short range potential charac- 
terized by the s-wave scattering length a. The 2D regime is 
realized when the trapping potential is sufficiently tight in one 
direction (which we take to be z) that Hu> Xty <C ksT ~ hut z , 
with T the system temperature. We limit our consideration to 
the regime where the scattering is three-dimensional (i.e. en- 
ergy independent), which requires that the z-confinement 
length, a z 



This require- 
11. It is conve- 



= \Jh/muj z , satisfies a z 3> a 
ments is well-satisfied by experiments (QUI 
nient to define two addition 2D regimes, which we will often 
refer to: The pure-2D regime, when ksT < ftio z , and the 
quasi-2D regime, when fcgT > fuv z . 



1. FSC theory 

If interactions are small compared to hw z then the Hartree- 
Fock modes are of the separable form ipka(x,y,z) = 
fka(x,y)^k{z), where the axial modes £k(z) are bare har- 
monic oscillator states. In the regime of interest the radial 
plane, for which we introduce the notation r = (x,y), can be 
treated semiclassically, eliminating the need to diagonalize for 
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the modes /fc CT (r) (also see lfl8tl ). However, the axial modes 
must be treated quantum mechanically, and the Hartree-Fock 
expression for the areal density of the system in the j-th axial 
mode is 
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V j( r ) = "r^tf + J huj z + 2 9kjn k (r), (3) 
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[i is the chemical potential, and 
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describes the interactions between atoms in the k and i axial 
modes. Performing the momentum integration in Eq. (Q]) and 
adding up the axial mode densities gives the total areal density 
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)=- ^m[l-exp {{n-V 3 (r)}/k B T)], (5) 
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where XdB = /i/V^rmfceT is the thermal de Broglie wave- 
length. We refer to the self-consistent solution of Eqs. (HJ and 
(O as the full semi-classical (FSC) mean-field theory, to dis- 
tinguish it from the simplified HCK theory presented below. 



2. HCK theory 



as where the central (peak) phase space density of the system 
satisfies the condition 
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with C = 380 ± 3 and g is the 2D coupling constant. Within 
the FSC theory the correct identifications for using Eq. © is 
g = goo and n cr = no(0) (i.e. the ground axial mode inter- 
action parameter and density, see discussion in 11611 ). however 
HCK proposed using g = ga and n cr = J2j n j (0) (i.e. total 
areal density). 

In combination, the differences in formulation of the two 
theories and their respective procedures for identifying Tbkt, 
is that the critical temperature predicted by FSC (Tj^p) is 
lower than HCK value (T^^). 



III. RESULTS 

To investigate these mean-field theories in the vicinity of 
the transition we compare their predictions to those obtained 
by a c-field method. The c-field method is based on a clas- 
sical field representation of the highly occupied low energy 
modes (i.e. those below an appropriately chosen energy cut- 
off, see M19I0 and a Hartree-Fock treatment of the sparsely 
occupied (high energy) modes of the system. This method 
treats the highly occupied (strongly fluctuating) modes non- 
perturbatively and is valid in the critical region. For details of 
this c-field theory refer to [20] and for the specific application 
to the quasi-2D trapped Bose gas see |0, 0, l2ltl . All results 
we present here are for 87 Rb atoms in cylindrically symmet- 



ric trap with radial trapping frequencies ui x 
Hz and axial frequency of either ui z 
uo z = 2tt x 1.88 kHz. 



j y — 2n x 9.4 
2tt x 0.94 kHz or 



The theory of HCK 1 1311 is a simplification of the FSC 
mean-field scheme presented above made by taking the 
interactions to be axial mode independent, i.e., simplify- 
ing the mean-field interaction term to 2 X^fclo 9fci n fc( r ) ~~ ¥ 
2gB_n(r), with the average interaction strength gn = 
(Airati 2 /to) J dzp(z) 2 , where p(z) is the density of a single 
atom in a harmonic oscillator of frequency ui z at temperature 
T. 



3. Comparison of theories and Tbkt 

An in-depth comparison the FSC and HCK mean-field the- 
ories is given in [16]. Most importantly, that study showed 
that the HCK theory predicts appreciably different density 
distributions in the degenerate regime, typically with higher 
(peak) densities. 

Mean-field theory can be used to estimate the superfluid 
transition temperature for the trapped system by using a local 
density extension of the transition criterion for the uniform 
system it 10TI . The transition temperature (Tbkt) is identified 



A. System density in the transition region 

An important quantity for detailed comparison is the sys- 
tem density profiles. In the development of the HCK theory 
direct comparisons of the density profiles were made to the 
results of quantum Monte Carlo calculations il3ll . however 
results were only given for cases where the peak phase space 
density was less than half the critical value ©, i.e. in rela- 
tively non-degenerate regime where the mean-field descrip- 
tion is expected to work well. Here we present results from 
regimes closer to the transition. 

Figure Q] shows density profile predictions of the various 
theories at several temperatures. The case shown in Fig.QJa) 
is approximately at the transition temperature (as determined 
by the c-field method [7]), and reveals a rather large difference 
between the c-field and mean-field predictions. In this regime 
both mean-field theories fail because they do not describe 
the quasi-condensate (suppression of density fluctuations) that 
forms [22]. Quasi-condensation allows the system density to 
greatly increase, explaining why the c-field density exceeds 
the mean-field density. The onset of quasi-condensation is 
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Figure 1: Phase space density. Parameters: (a) , N = 13.4 x 10 3 , 
T = 34.7 nK (« Tbkt). (b) N = 45.8 x 10 3 , T = 67.6 nK. (c) 
N = 67.4 x 10 3 , T = 65.0 nK (« Tbkt). Curves: FSC (solid), 
HCK (dashed), and c-field (dot-dashed) results. Other parameters: 
Lo x , y /2-iT = 9.4 Hz for all results, w z /2tt = 1880 Hz in (a) and (b), 
and u) z /2n = 940 Hz in (c). 



not abrupt at the transition, but gradually forms as a precur- 
sor to the transition at temperatures well above Tbkt ifTlll 
(e.g. in Fig. [T}b) where a small quasi condensate component 
[about 3.6% of the system, but quite dominate near trap cen- 
ter] causes the c-field density to be slightly higher than the 
mean-field predictions near r = 0). The results in Fig. [Tfa) 
and (b) are in the pure-2D regime, and in Fig.[TJc) we consider 
a case near Tbkt, but in the quasi-2D regime. For this case, 
the appreciable thermal activation in the z-direction reduces 
<7h from 500 an d the HCK theory is closer to the c-field result. 



B. Comparison of predictions for Tbkt 

We can now compare the predictions of the FSC and 
HCK mean-field theories for the transition temperature Tbkt 
against the c-field simulations. In Fig. [2] we present an exam- 
ple of how this comparison is made. For this case the system 
has a large aspect ratio (i.e. uj z /uj x . y — 200) and in the transi- 
tion region the system is in the pure 2D limit where goo ~ 3h, 
and the difference between no (r) and the total areal density is 
negligible. 

For the two mean-field theories we calculate the peak phase 
space density at each temperature and then locate the tran- 




Figure 2: Peak phase space density as a function of T. FSC value 
of no(0)A^ s (pluses), HCK peak value of the total areal phase space 
density, ]TV n j W^dB (triangles) and polynomial fits to these results 
in the transition region are shown as lines. Also displayed (as nearly 
horizontal solid lines) are the critical phase space densities for the 
FSC and HCK theories (see Eq. (JSJl). The predictions for the transi- 
tion temperature made by the FSC {T£ K §), the HCK (T§gfi) and 
c-field (Tbkt) theories are indicated as vertical dashed lines. Cal- 
culation parameters: u x ,y = 2ir x 9.4 Hz, u z ,= 2n x 1.88 kHz and 
N = 13.5 x 10 3 . 



sition temperature by condition © appropriately interpreted 
for each theory (see discussion in Sec. Ill 31 . The respective 
critical phase space density for each theory is shown Fig. [2] 
however for the pure-2D regime of this calculation there is lit- 
tle difference between these values, and we find that Tf|g 
and TggT are similar to each other but are both significantly 
reduced from the c-field prediction for Tbkt- 

In Table U we repeat the analysis shown in Fig. [2] for a 
system in five different parameter regimes, which we refer 
to as cases A-E. Importantly, these results explore parame- 
ters whereby the system crosses over from being pure-2D to 
quasi-2D at the transition point, as quantified by the parameter 
v z = ksTsKT /foj->z- In the pure- 2D regime (v z < 1, cases A, 
B) we find that the FSC and HCK theories make similar pre- 
dictions, and both significantly underestimate the transition 
temperature. This occurs because they fail to account for the 
suppression of density fluctuations and thus underestimate the 
phase space density of the system, as observed in Fig.[Tfa). In 
the quasi-2D regime (v z > 1, cases C, D, E) the HCK predic- 
tion agrees quite well with the c-field result for Tbkt, while 
the FSC remains an (appreciable) underestimate. The main 
origin of the improved agreement is the reduction in gjj that 
occurs as v z increases, which tends to increase the HCK den- 
sity. We emphasize this increase is not for the correct physical 
reason (suppressed density fluctuations), but simply from the 
ad hoc choice of gjj- Figure [T](c) shows the density profile 
near Tbkt for case D. We also note that Tbkt is identi- 
fied in the c-field simulations using the same criterion used 
for the FSC theory, i.e. when the peak phase space density of 
the ground axial mode satisfies the condition (O with g = goo- 
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Cases 




A 


B 


C 


D 


E 


AT/10 3 


13.5 


21.2 


30.2 


67.3 


228.5 


u> z /2tt 


1880 Hz 


1880 Hz 


940 Hz 


940 Hz 


940 Hz 


Tbkt 


34.7 nK 


43.1 nK 


46.8 nK 


64.9 nK 


104 nK 




0.38 


0.48 


1.04 


1.44 


2.32 


rpFSC 

1 BKT 


31.0 nK 


38.4 nK 


43.5 nK 


60.4 nK 


96.9 nK 


rel. err. 


n i 1 
-U.l 1 


n 1 1 

-U. 1 1 


-U.U / 


-U.U / 


-U.U / 


m ( 380ft 2 \ 


8.18 


8.18 


8.52 


8.52 


8.52 


rpHCK 

1 BKT 


31.4 nK 


39.2 nK 


45.8 nK 


64.4 nK 


105 nK 


rel. err. 


-0.10 


-0.09 


-0.02 


-0.01 


+0.01 


In ( 380h2 ) 

V m SH / 


8.23 


8.28 


8.92 


9.07 


9.30 



Table I: Comparison of mean-field predictions for the critical tem- 
perature. The c-field prediction for the critical temperature {Tbkt) 
is used to obtain v z . The relative error of each mean-field estimate 
of Tbkt is given against the c-field value and the radial trap fre- 
quencies in all cases are ui XtV /2n — 9.4 Hz. Interaction strength 
qh calculated using corresponding temperature Tg^T ■ Logarithms 
indicate critical phase space density condition for each theory. 

In Ref. lfl3ll the comparison of the HCK prediction for 
Tbkt with QMC is made for a system with v z « 1.4, and 
agreement is found to be at the few percent level. In the con- 
text of our results here, it is clear that this good agreement 
is likely facilitated by the reduction in gjj from 500 compen- 
sating for the absence of quasi-condensate in the mean-field 
theory. 

IV. CONCLUSIONS 

We have compared the density predictions made by two 
semiclassical mean-field theories against the c-field method. 



At low temperatures, we find that both mean-field theories 
overestimate interaction effects due to their neglect of cor- 
relations that cause density fluctuation suppression (quasi- 
condensation). In this regime both methods systematically 
underestimate the central density, and underestimate Tbkt- 

For larger systems, where the transition occurs at tempera- 
ture with appreciable thermal excitation in the tightly trapped 
mode, the scenario is quite different. The HCK approximation 
of replacing the mode dependent interaction parameter by a 
thermally averaged value, leads to an increase in the system 
density, and a higher prediction for the transition temperature 
than the FSC theory. We expect that for sufficiently large sys- 
tems the HCK theory will over estimate Tbkt, although we 
have not presented results here that show this. 

The quantitative tests we have presented provide a useful 
characterization of the applicability of mean-field methods in 
the vicinity of the BKT transition for the 2D trapped Bose 
gas. Our results are over a broader range of parameters than 
those considered in O, and thus present a more balanced 
view of the method effectiveness. Generally, our results show 
that mean-field methods cannot be relied upon for accurate 
predictions, but do nevertheless provide a rough estimate of 
the degeneracy regime. With future experiments likely pro- 
ducing tighter 2D traps, in which degeneracy will occur in the 
pure-2D regime, our observations on the failure of mean-field 
theory in this regime will be important. Of course, in the quest 
for a simple and effective mean-field theory several other fac- 
tors must be borne in mind: (i) The HCK theory involves a 
single interaction parameter allowing it to be computed ex- 
tremely efficiently using the procedure presented in II 1311 . (ii) 
The FSC theory provides a lower bound for Tbkt- 
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